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ABSTRACT 

Classical novae are the result of thermonuclear flashes of H accreted by CO 
or ONe white dwarfs, leading eventually to the dynamic ejection of the surface 
layers. These are observationally known to be enriched in heavy elements, such 
as C, O and Ne that must originate in layers below the H-flash convection zone. 
Building on our previous work we now present stellar evolution simulations of 
ONe nova, and provide a comprehensive comparison of our models with previ- 
ous work. Some of our models include exponential convective boundary mixing 
model to account for the observed enrichment of the ejecta even when accret- 
ing material with a solar abundance distribution. Our models produce maximum 
temperature evolution proflles and nucleosynthesis yields in good agreement with 
models that generate enriched ejecta by assuming that the accreted material was 
pre-enriched. We conflrm for ONe nova the result we reported previously, i.e. 
we found that ^He can be produced in situ in solar-composition envelopes ac- 
creted with slow rates (M < 10-^° Mq/jt) by cold (Two < 10^ K) CO WDs, 
and that convection is triggered by '^He burning before the nova outburst in this 
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case. In addition, we now find that the interplay between the ^He production 
and destruction in the solar-composition envelope accreted with an intermediate 
rate, e.g. M = 10'^^ Mq/jt, by the 1.15 Mq ONe WD with a relatively high 
initial central temperature, e.g. Two = 15 x 10^ K, leads to the formation of a 
thick radiative buffer zone that separates the bottom of the convective envelope 
from the WD surface. We present detailed nucleosynthesis calculations based 
on the post-processing technique, and demonstrate in which way much simpler 
single-zone trajectories extracted from the multi-zone stellar evolution simula- 
tions can be used, in lieu of full multi-zone simulations, to analyse the sensitivity 
of nova abundance predictions on nuclear reaction rate uncertainties. Trajecto- 
ries for both CO and ONe for different central T and accretion rates are provided. 
We compare our nova simulations with observations of nova and pre-solar grain 
believed to originate in nova. 

Subject headings: methods: numerical — stars: novae — stars: abundances — 
stars: evolution — stars: interiors 



Introduction 



Classical nova explosions are a consequence of thermonuclear runaways (TNR) occurring 
in accreted H-rich shells on the white dwarf (WD) components of close binary systems. 
The accumulation of H-rich shells on the surfaces of the WD components of these systems 
continues until a critical pressure ~ lO^^dyne/cm^ is achieved at the base of the accreted 
envelope and runaway ensues. An increase in the shell burning luminosity then follows (at 
a rate that reflects/defines the speed class of the nova) to a value that approaches (and for 
fast novae can exceed) the Eddington limit. Following a period of burning at app roximately 



constant bolometric luminosity, dictated by the core- mass-luminosity relation of iPaczyhski 



(jl970l ). the nova ultimately returns to its pre-outburst state. The occurrence of mass loss 
from novae through outburst yields an ejected envelope, the composition of which reflects 
both the consequences of the thermonuclear burning epoch and the occurrence of dredge-up 
of matter from the underlying WD into the H-rich envelope. 

Our understanding of the nova phenomenon has increased substantially over the past 
several decades as a consequence of significant progress in theoretical modeling of their out- 
bursts, of significantly improved observational determinations of the characteristics of novae 
throughout outburst, and of the acquisition of greatly improved abundance data concerning 
both the comp ositions of th e accreted shells and, particu l arly, the nebular e jecta of diverse 
nova systems (ITruraru Il982l : IStarrfieldl Il989l : ISharal Il989l : iGehrz et al.lll998l ). Multidimen- 
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siona l hydro dynamic simulations of classical novae are also available (iGlasner et al.l Il997 
20051 . l2012[ ICasanova et~aDl201o[ boilbl : iKercek et allllQQsh . 



A critical feature of novae seen in outburst is the fact that diverse nova systems have 
been found to be characterized by significant enrichments of many of the elements, such as 
C, N, O, Ne, Na, Mg, or Al, relative to H, with respect to solar abundances. It is generally 
expected that, for the temperature/density conditions that characterize the burning shells of 
novae (with peak shell temperatures typically less than 4 x lO^K), significant breakout of the 
CNO H-burning sequences cannot occur. It then follows that the elevated abundance levels 
arise from dredge-up from the underlying cores. C, N, and O are enriched by dredge-up from 
CO WDs, while the observed abundance enrichments of C, O, Ne and maybe Mg arise from 
the outward mixing or dredge-up of matter from an underlying ONe WD. This picture is 
to be contrasted by the commonly adopted practise that the accreted envelope in numerical 
studies is pre-enriched in order to mat ch the observation s of the abundances in nova nebular 
remnants (which are reported in, e.g. iGehrz et al.lll998l ). 



ONe nova are of particular interest for nucleosynthesis because they develop the highest 
temperatures during the nova outburst. They also represent an unexpectedly large fraction 



of observed novae. Tabl es of abundances in no va ej ect a were publis 



on nova composition b y iTruran fc Livid (119861 ) and iLivio fc Truran 



l ed in both the papers 
( 1994 ) and the review 



paper by lGehrz et al.l (jl998l ). In the more extended compilation by iGehrz et al.l ( 1l998l ) out 
of the 20 novae ten are clearly CO novae (Ne/Ne© < 8) while seven (three) are clearly 
(possibly) ONe novae (Ne/Ne0 > 10 — 15). However, the fraction of ONe white dwarfs 
in the Galaxy is expected to be much smaller, based on models of stellar evolution. C 
ignition happens only for (super-)Asymptotic Giant Branch (AGB) core masses greater than 
~ 1.05Mq, and ONe WD can only form from stars with initial masses in the corresponding 
narrow initial mass range ~ 7 — IOMq (depending on metallicity and mixing assumptions 
during H- and He-core burning phases) that features C ignition but avoids core-collapse 
supernova. Although the details of this sim ulation prediction depends sensitively on the 
i2q _|_i2 q reaction rate (IChen et al.ll2013() as well a s the physics details of the C-fiame 
for off-center C ignition (IDenissenkov fc Herwig] l2013l ). present best estimates based on a 
Salpet er intial mass function predict that the ra tio of CO to ONe WD is approximately 
thirty Jhivio fc Truranlll994l : Icil-Pons et al.ll2003[ ). 



Denissenkov et al.l (120121 . hereafter Paper I) presented ID simulatio ns of nova outbursts 



occurring on CO WDs based on the MESaEI stellar evolution code (iPaxton et al.l l201ll 



20131 ). In the present paper, we show present MESA models of the ONe novae outbursts. In 



^MESA is an open source code; 



http : //mesa. sourceforge .net 
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addition, th e multi-zone post-pr ocessing; code MPPNP, which is a part of the NuGrid research 



framework (IHerwig et al.ll2008l ). is used here for computations of detailed nucleosynthesis 
in CO and ONe novae. The MESA nova models provide MPPNP with necessary stellar 
structure and mixing parameter profiles. The NuGrid tools also include the single- zone code 
SPPN that enables post-processing nucleosynthesis computations along a (T, p) trajectory, 
for example at one Lagrangian coordinate (a mass zone) inside ID stellar models, or along 
the position of the maximum temperature. We have run SPPN for trajectories tracing the 
evolution of the maximum temperature in the nova envelope (usually at the core-envelope 
interface) and its corresponding density. We show that, because of the strong dependence 
of thermonuclear reaction rates on T, nucleosynthesis yields obtained with the SPPN code 
for the Tmax-trajectories are in qualitative agreement with the complete yields provided by 
the MPPNP code. The SPPN code runs much faster than MPPNP, therefore it can be used 
for a comprehensive numerical analysis of parameter space, for example, when the relevant 
reaction rates are varied within their observationally and experimentally constrained limits. 

We have prepared shell scripts, template MESA inlist files, and a large number of initial 
WD models aiming to combine MESA and NuGrid into an easy-to-use Nova Framework. 
This new research tool can model nova outbursts and nucleosynthesis occurring on CO and 
ONe WDs using up-to-date input physics and a specified nuclear network for a relatively 
large number (up to a few thousands) of mass zones covering both the WD and its accreted 
envelope. Nova simulation results can be analyzed using animations and a variety of plots 
produced with NuGrid's Python visualization scripts. One of the main goals of this paper 
is to demonstrate Nova Framework's capabilities. 

The simulation tools of MESA and NuGrid are described in Section [2l The MESA mod- 
els of nova outbursts occurring on CO and ONe WDs are discussed in Section [31 Results 
of the post-processing nucleosynthesis computations for CO and ONe novae obtained with 
the NuGrid codes MPPNP and SPPN are presented in Section HI In Section [S], we describe 
the effects of ^He in situ production and burning on triggering convection in the accreted 
solar-composition envelopes of ONe novae. A comparison of the results of the nova nucle- 
osynthesis simulations with available data on nova-related element and isotope abundances 
from spectroscopic and presolar grain composition analyses is made in Section O Finally, 
our main conclusions are drawn in Section [71 
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Simulation Tools 



2.1. Stellar Evolution Calculations - MESA 



MESA is a collection of Fortran-95 Modules for Experiments in Stellar Astrophysics. 
Its star module can be used for ID stellar evo l ution s imula tions of almost any kind (nu- 
merous examples are provided by iPaxton et al.l (120111 . |2013| ) and on the MESA website). 
MESA modules provide star with state-of-the-art numerical algorithms, e.g. for adaptive 
mesh refinement and timestep control, atmospheric boundary conditions, and modern input 
physics (opacities, equation of state (EOS), nuclear reaction rates, etc.) 

For our MESA simulations of ONe nova outbursts, we use the same EOS and opacity 
tables that were used to model the CO novae in Paper 1. The only important difference is 
our adoption of longer lists of isotopes and reactions for the ONe nova models. The adopted 
nuclear network has to be as large as necessary in order to account for the energy generation 
during a nova TNR, yet as small as possible in order to make computations not too expensive. 
To make an optimal choice, we have started our MESA nova simulations with the largest 
network nova_weiss .net that includes 77 isotopes from H to ^'^Ca coupled by 442 reactions. 
This n etwork has a li s t of i sotop es almost identi c al to that selected by IWeiss fc Truran 
(Il990l ) , iPolitano et al.l (119951 ) , and IStarrfield et al.l (120091 ) . We have gradually reduced the 
numbers of isotopes and reactions checking that this does not lead to significant changes in 
the Tinax-trajectories. As a result, acceptable compromises have been found empirically. A 
list of 33 isotopes coupled by 65 reactio ns, including tho s e of t he pp chains (the pep reaction, 
whose importance was emphasized by IStarrfield et al.l (120091 ) . has also been added), CNO 
and NeNa cycles, was selected for our CO nova models in Paper I. We will call this CO 
nova nuclear network simply nova.net. A similar procedure for the ONe nova models has 
resulted in a necessity to adopt an extended nuclear network nova_ext . net that includes 
48 isotopes from H to '^°Si coupled by 120 reactions. The Tmax-trajectories extracted from 
our CO and ONe nova simulations performed for different sets of nova parameters using 
the three nuclear networks are compared in the upper panels of Figs.iH IH and [HI From 
these plots, it is seen that the nova . net and nova_ext . net nuclear networks represent an 
acceptable accuracy for the simulations of CO and ONe nova outbursts, respectively. The 
largest network nova_weiss .net can optionally be used for extreme cases and test purposes. 



By default, MESA uses reaction rates from lCaughlan fc Fowlerl (119881 ) and lAngulo et al. 



(119991 ) ■ with preference given to the second source ( NACRE). It includes up dates to the 
NACRE rates for ^^ N(p,7)^^0 (jlmbriani et al.l |2005| ) , the t riple-a reaction (iFvnbo et al. 



20051), ^^N(a,7)iSF dCorres et al.l l2000h . and ^^C{a,'yy'^0 dKunz et al.ll2002h . Although 
the main nuclear path for a classical nova is driven by p-capture reactions and /3-decays, 
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the a-reactions are important for establishing the chemical composition of its underlying 
WD. As a test, we have also tried t he MESA option fo r choosing the reaction rates that 
gives preference to the JIN A Reaclib (ICvburt et al. I20101). This datab ase has many nuclear 
reaction rates in common with those considered by llliadis et al.l (120101 ) . We have not found 
any significant differences in our nova Tmax-trajectories between the two options (e.g., see 
the upper panel in Fig. |6]), therefore we decided to stick to the default one. 



2.2. Nucleosynthesis Calculations - NuGrid 



The nucleosynthesis calculations were performe d with the post-processing code PPN (ei- 
ther for single zones or multi-zone post-processing, iHerwig et al.ll2008t iPignatari fc Herwig 
20121 ) . where the input data of stellar structure (here from MESA models) are processed 
using a dynamically updated list of nuclear species and reaction rates. The network can 
include more than 5000 species, between H and Bi, and more than 50000 nuclear reactions. 
The self-adjusting dynamical network has at any point in the simulation (and in each layer 
for multi-zone simulations) the optimal network size, based on the strength of nucleosynthe- 
sis fluxes. Rates are tak en from different da ta sources: the European NACRE compilation 



2005 



(lAngulo et al.lll999[) andllhadis et al.l (20011) . or more recent, if available fe.g. . iFvnbo et al. 



Kunz et al.l l2002l : llmbriani et al.l 120051 ) , and the JINA Reaclib library (ICyburt et al. 



20101 ). For neutron capture rates by stable nuclides and a selection of relevant unstable iso- 
topes, we use the Kadonis comp i lation (http : / / www . kadonis . org) . For weak rea c tion r ates, 
we choo s e betw een iFuUerl ( 119851 ) , lOda et al.l (Il994j ) , iLanganke fc Martfnez-Pinedol (120001 ) and 



Gorielyl (119991 ). according to the mass region. 



For the complete post-processing of the full MESA nova models, we use the multi-zone 
parallel frame of PPN (MPPNP). Trajectories of single (T, p) evolutions are computed with 
the single-zone version SPPN. Both variants use the same nuclear physics library and the 
same package to solve the nucleosynthesis equations. 
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3. MESA Models of Nova Outbursts 

In Paper I, we presented CO nova outburst MESA models for a number of initial CO 
WD masses and temperatures (luminosities). For completeness we include CO nova burst 
models here as well and then present the new ONe nova models. The nucleosynthesis in 
both model families is then analysed using the NuGrid tools. 



3.1. CO Novae 

A few CO nova models are included in the present discussion as a link to Paper 1 and also 
as an illustration of application of the NuGrid MPPNP code for post-processing of CO nova 
models. All the results relevant to the CO nova models are presented in Figs. [T]and[T21 and at 
the beginning of Tabled! The upper panel in Fig. [T] compares two Tmax-trajectories extracted 
from our MESA simulations of nova outbursts occurring on a CO WD with the same mass, 
MwD = 1.15 M0, and initial central temperature, Two = 10 MK (log^Q Lwd/-Z^0 = —2.69; 
for the correspondence between WD's initial central temperature and luminosity. Two and 
LwD; see Tabled]), and for the same accretion rate, M = 10~^^ Mq/jt. The only difference 
between the models is that the one has been computed using the small nuclear network 
nova.net that includes 33 isotopes from H to ^^Mg coupled by 65 reactions, while the 
other with the extended network of 48 isotopes and 120 reactions, nova_ext .net. The good 
agreement between the two trajectories, even for a relatively low Two and slow accretion 
rate that leads to a stronger outburst, confirms our choice of nova.net as a sufficient nuclear 
network for the MESA CO nova models in Paper 1. The lower panel shows the results of 
the post-processing nucleosynthesis simulations carried out for one o f the CO nova rnodels 



considered in Paper I. Our results qualitatively agree with those of (iJose fc Hernang 11998 
see their Fig. 1) for a CO nova model with similar parameters and with those produced by 
the MESA code alone (see Fig. 5 in Paper I). 



3.2. ONe Novae 



According to iGil-Pons et al.l (120031 ). the frequency of occurrence of galactic ONe novae 
in close binary systems is expected to be 30% to 40%, depending on model assumptions. 
Without convective overshooting of H-core burning, the range of initial masses of the stars 
that end their evolution as ONe WDs with the masses between I.IMq and 1.33Mq is esti- 
mated by the same authors to be between 9.3M0 and 12M0. 



Two ONe WD models with the masses Mwd = 1.15 Mq and 1.3 Mq have been prepared 
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for a range of initial central temperatures and luminosities using the same "stellar engineer- 
ing" procedure with which we created the CO WD models in Paper I (Table [1], note that 
the Lwd(^wd) relation is almost universal for the CO and ONe WDs for Mwd > I.OMq). 
Likewise, we have artificially removed the He- and C-rich buffer zones from the surfaces of 
ONe WD models and used the naked ONe WDs with thin (~ 10"'' Mq) H-rich envelopes 
as the initial models for ONe nova simulations. The ONe WD models were created using 
the same isotopes as in nova_weiss.net, while the reaction list was extended to take into 
account He and C burning. 

We have not included any prescription for mass loss in the MESA ONe nova simulations, 
yet we have followed the hydrodynamic expansion of the nova envelopes until their surface 
radii exceeded a few solar radii. This should guarantee that the expanding envelope has 
crossed the WD's Roche- lobe radius for a typical binary system hosting a classical nova 
(e.g.. Paper I). The ONe nova evolutionary tracks are not much different from our CO nova 
tracks discussed in Paper I (e.g.. Fig. [2]). A more detailed discussion of nova model light 
curves is beyond the scope of the present study, which is mainly focused on the mechanism 
of the outburst phase and the nova nucleosynthesis. 

Like for the CO nova simulations in Paper I, here we consider three types of ONe nova 
models. A few models are of the first type and neglect boundary mixing (adopting a strict 
Schwarzschild condition to determine the convective boundary) and assume that the WD ac- 
cretes solar-composition material (e.g., the models with the in situ ^He production presented 
in Figs. [T71 and [T8|) . These models do neither show the fast rise time nor the enhancements 
in heavy elements observed for nova. The second type take convective boundary mixing 
(CBM) (e.g., the CBM models in Figs.fTT] and [T2|) into account, via a diffusion coefficient 
that is exponentially decreasing with geometric distance below the bottom of the convective 
envelope. The corresponding e-folding length is / ■ Hp, where Hp is the pressure scale height 
near the Schwarzschild boundary, and / = /nova = 0.004 has the same value as in Paper I. 
These models do match observed rise times and abundances. 

Finally, most of our models (Table [1]) are of type three in which CBM is neglected, like 
in the first-type models. However, the effect of CBM is artificially obtained by assuming that 
the WD accretes a mixture in which solar-composition material has already been blended 
with an equal amount of WD material. 

Such pre-mixed models are commonly adopted in ID nova simulations. The fraction 
of W D's material in the m ixture can vary from 25% up to 75% to match observed values 



le.g- 



Jose fc Hernandll998l ). To prepare the initial pre-mixed abundances, we have used the 
isotope abundances from the outermost layers of our naked ONe WD models. Our initial 
pre-mixed abundances for the I.SMq ONe nova model (of the third type) are compared with 
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the 50% pre-mixed abundances that we c all the Bar c elona and P olitano compositions in 
Fig. [3l The Barcelona data are taken from Jose et al. ( 1999 . 200 ih. They contain m aterial 
from a l.SSM© ONe WD model. The Politano data are from iPolitano et al.l (11995 ). The y 
were derived from the C-burning nucleosynthesis calculations of lArnett fc TruranI ( 1l969l ). 
Our initial abundances agree quite well with those from the Barcelona composition (the 
upper panel), especially for ^^O, ^°Ne, and ^^Mg. As for the Politano composition (the lower 
panel of Fig. [3]), it has similar to ours ^^O and ^"Ne abundances, but much higher initial 
^^Mg abundance. The last difference is important for ONe nova simulations because ^^Mg 
can compete with ^^C in ign iting the accreted H -rich material in the range of temperatures 
relevant to nova outbursts flGlasner et al.l |2012| ). Our initial ^^C abundance has a value 
intermediate between the two others. 

In Fig.m we compare Tmax-trajectories extracted from our I.ISMq (panels A and B) and 
1.3Mq (panels C and D) ONe nova simulations for Two = 12 MK with pre-mixed material, in 
which isotope abundances from WD models of the corresponding masses were used, accreted 
with the rate M = 2 x 10-^°Mo/yr. For either CO or ONe WD, the main properties 
of its nova outburst, such as total accreted (Mace) and ejected masses, peak temperature 
Tjnax, maximum H-burning (Lh) and total luminosities, envelope expansion velocity, and 
chemical composition of the ejecta, depend mainly on the following three parameters: the 
WD i nass Mwp, its central ternperature Twn (or luminosity Lwn), and the accretion rate M 



(e.g., iPrialnik fc Kovetzlll995l : iTownsley fc Bildstenll2004l : I Jose et al.ll2007l ). The results also 
depend on a choice of nuclear network (lists of isotopes and reactions), initial abundances in 
the accreted material, reaction rates and other input physics (EOS, opacities), and computer 
code. For the present work, some of the initial and calculated nova model parameters 
are listed in Table [1] Here, we have only considered material (and mixtures thereof) of 
solar metallicity. In general, our pre-mixed ONe nova models computed with the Barcelona 
initial composition have systematically lower M^^c and Tmax values, as compared to their 
counterparts computed with both MESA and Politano compositions. This is probably caused 
by the higher initial abundance of ^^C in the Barcelona mixture (the upper panel of Fig. [3]) 
that turns on the nova ignition a bit earlier. 

To categorize our nova simulations carried out for the different sets of initial parameters, 
we use Tjnax-ti'ajectories (e.g.. Fig. HJ. They have two characteristics that determine the 
resulting nova speed class and the nucleosynthesis: a peak value of the Tmax curve and its 
width near the peak, the latter being proportional to a timescale of nova outburst. We 
have already used Tmax-trajectories to support our choice of nova.net as a sufficient nuclear 
network for MESA CO nova simulations (Fig. (H upper panel). Similarly, panels A and B of 
Fig. m show that our extended nuclear network is acceptable for MESA ONe nova simulations 
because its further extension to nova_weiss.net does not change the trajectory (panel B), 
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while using of nova.net results in a slower nova outburst (panel A). Panel C demonstrates 
that changing the initial composition in the accreted envelope from the Barcelona to Politano 
one produces differences in the Tmax-trajectories comparable to those caused by the replacing 
of nova.net by nova_ext .net in panel A. The small difference between the curves in panel 
D confirms the conclusion derived from Fig. |3]that our pre-mixed initial composition is closer 
to the Barcelona one. 

The Tjnax-trajectories and their corresponding density profiles can be used by the NuGrid 
SPPN code for post-processing computations that trace nucleosynthesis for the specified 
variations of T and p with time. A shortcoming of such computations is that they do not 
take into account abundance changes occurring in other parts of the nova envelope that 
are connected by convective mixing. Much more consistent post-processing nucleosynthesis 
simulations can be done with the NuGrid MPPNP code that uses detailed radial profiles 
of T, p and the diffusion coefficient from the stellar evolution simulation to account in the 
same way for time-dependent mixing. Our Nova Framework uses a MESA option that allows 
to generate such input stellar models for MPPNP, during computations of nova evolution. 
Still, given that thermonuclear reaction rates strongly depend on T, even one-zone SPPN 
simulations should probably produce results that are, at least qualitatively, compatible with 
those obtained in full MPPNP simulations, as discussed in more detail in the next two 
sections. 



4. Post-Processing Computations of Nova Nucleosyntheis 

4.1. Multi-Zone Computations 

MESA simulations are using a limited nuclear network. The limit is imposed by the 
number of isotopes and reactions needed to accurately simulate the nuclear energy gener- 
ation rate which, at the same time, should not be too large to make MESA computations 
prohibitively slow. The NuGrid multi-zone code MPPNP allows to study nucleosynthesis 
in stars more accurately than MESA because it can include nuclear networks with more 
than 5000 isotopes coupled by more than 50000 reactions (see § 12. 2^ and references therein). 
This is possible because MPPNP solves for mixing and burning in an operator split while 
MESA, which also has to solve for the structure and energy equations benefits from the 
much improved convergence properties of a fully coupled operator approach. MPPNP works 
in a post-processing mode using stellar structure models (cycles) prepared by MESA as a 
input. For nova nucleosynthesis simulations, we have identified a list of 147 isotopes coupled 
by more than 1700 reactions. This nuclear network includes all the species relevant for nova 
nucleosynthesis (Fig. [5]). As an example, we show the magnitudes of the main reaction fiuxes 
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for the trajectory with the largest temperature (Tmax = 408 MK) in Fig. |5l 

In Fig. El upper panel, we show the Tmax-trajectories from two I.SMq ONe nova models 
with = 20 MK that accreted the pre-mixed material with the rate M = 2 x 10~^'^MQ/yr 
for different choices of reaction rate sources available in MESA, NACRE and JINA Reaclib. 
The two trajectories do not differ much from one another (see the discussion in the previous 
section). In particular, their corresponding nuclear energy generation budget differs only by 
16%. 

The lower panel of Fig. [6] shows the production factors (mass fractions normalized to 
the solar composition) for stable isotopes whose abundances were averaged over the mass of 
the expanding nova envelope for the two discussed nova models. The nucleosynthesis yields 
are comparable in the two cases, with the largest differences for ^^Ne and Mg isotopes. In 
the lower panel of Fig. [H the results of our MPPNP simulations are shown for one of the 
CO nova models from Paper I for which a similar plot, but with data extracted from MESA 
simulations, was displayed there in Fig. 5. The nucleosynthesis yields are presented in the 
same form and for a CO nova model with the similar parameters in Fig. 1 of I Jose fc Hernanz 



( 1l998l ) as well. The data in the three figures demonstrate a very good qualitative agreement, 
especially given that they were obtained with different nuclear networks, reaction rates, 
initial isotope abundances, and stellar codes. 

The final abundances from four nova models are shown in Figs. [7] and |H1 Two models 
are characterized by high WD temperature and relatively fast accretion, calculated using 
two different networks (nova_wiess.net and nova_ext.net. Fig. [7]). Their Tmax-trajectories 
and final abundances are really similar, confirming the result in Fig. HI panel B. For the 
nova models with the coldest WDs and accreting with the slowest rates (Fig. [8]), the Tmax- 
trajectories are more different, although the relative nova yields do not show significant 
variations between the two cases. In Fig. [HI we plot our final MPPNP abundances (upper 
panel) in the I.SMq ONe nova model with Two = 12 MK and M = 2 x IQ-^^Mq/ji. In 
this case, both the MESA and MPPNP simulations have used the Politano composition 
for the accreted material (the lower panel of Fig. |3]). In the lower panel of Fig. |9l we 
show the ratios of our abundances to t he correspondin g final abundances from the ONe 



nova model 12005 that was computed by IStarrfield et al.l (120091 ) for the similar parameters 



in cluding the same pre-m ixed initial composition. Most of the abundances agree with those 



of IStarrfield et al.l ( l2009l ) withi n a factor of ten . A si milar comparison for our ONe nova 



model and its counterpart from iJose fc Hernanzl (119981 ) in Fig. [TOl both computed with the 



Barcelona initial composition (the upper panel of Fig. EJ, leads to the same conclusion. The 
final abundances from the two independent works with which we compare our results differ 
among each other as much as each differs from our results. A more detailed comparison 
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would require simultaneous access to all participating simulation tools, and may not be 
justified given the present accuracy of available observable constraints. 

A realistic nova model should assume that the accreted material has solar composi- 
tion (for solar metallicity) and that it is mixed with the WD's material before or during 
its thermonuclear runaway. Recent two- and three-dimensional nuclear-hydrodynamic sim- 
ulations of a nova outburst have shown that a possible mechanism of this mixing are the 
hydrodynamic instabilities and shear-flow turbulence induced by steep ho rizontal velocity 



gradients at the bottom of the convection zone triggered by the runaway (jCasanova et al. 



2010l . l2011al Jbl). These hydrodynamic processes associated with the convective boundary lead 
to convective boundary mixing (CBM) at the base of the accreted envelope into the outer 
layers of the WD. As a result, CO-rich (or ONe-rich) material is dredged- up during the 
runaway. In our ID MESA simulations of nova outbursts, we use a simple CBM model that 
treats the time-dependent mixing as a diffusion process. The model approximates the rate 
of mixing by an exponentially decreasing function of a distance from the formal convective 



boundary (IFreytag et al.lll996l : iHerwig et al.lll997l ) 



(2 1 ^ — I \ 
]Hp ) ' 

where Hp is the pressure scale height, and Dq is a diffusion coefficient, calculated using a 
mixing-length theory, that describes convective mixing at the radius ro close to the boundary. 
In this model / is a free parameter, for which we use the same value / = /nova = 0.004 that 
we used in our CO nova simulations in Paper I. For a physical motivation to model the 
CBM with the prescription ([T]), see Paper I. Our MESA nova models with the CBM have 
also been post-processed with the MPPNP code, and the resulting final abundances are 
compared with those obtained for the corresponding nova models without CBM but with 
the pre-mixed accreted envelopes in the lower panels of Figs. [11] and [12] for the ONe and 
CO novae, respectively. The comparison between the final abundances shows a very good 
agreement for both the ONe and CO nova models. This gives a support to the widely-adopted 
ID nova model in which the CBM is mimicked by assuming that the accreted envelope has 
been pre-mixed with the WD's material. 

In Fig. [131 we make a comparison between three models where the WD initial central 
temperature (luminosity) is varied. Significant differences are obtained over all the nova 
yield distribution. In particular, the final abundances increase with increasing the central 
temperature. 
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4.2. Single-Zone Computations: a Tool for Uncertainty and Sensitivity Studies 

In this section, we discuss the post-processing simulations that use Tmax-trajectories 
extracted from the MESA I.3M0 pre-mixed ONe nova models (shown in Fig. [T3|) with the 
mass accretion rate M = 2 x 10~^*^MQ/yr, but for the three different WD initial central 
temperatures. Two = 12, 15, and 20 MK . The accreted material had the Barcelona initial 
composition. The Tmax peak values for these models are 344, 313, and 267 MK, respectively 
(Table [1]). The post-processing nucleosynthesis calculations for these trajectories are per- 
formed using the NuGrid SPPN code. The solver and nuclear physics packages are the same 
that were used in our multi-zone post-processing nova computations ( §4.ip . and we use the 
same nuclear network with 147 isotopes (Fig. The final abundance ratios for stable iso- 
topes are presented in Fig. [TH The abundances were scaled using their corresponding solar 
values (the lower panel) as well as the results obtained for the WD model with Two = 20 
MK (the upper panel). 

A comparison between the abundance distributions in Figs. [T^ and [TUl shows a qualitative 
agreement between the results from the multi-zone complete models and from the single- 
zone Tmax-trajectories. Indeed, the abundance patterns for different groups of isotopes look 
similar. The main quantitative difference is the more efficient production toward heavier 
species in the second case, where O is more depleted, while Ar and K isotope^ are more 
efficiently made. This difference is caused by the fact that in the multi-zone post-processing 
simulations convective mixing reduces an "average temperature" at which the nucleosynthesis 
occurs and it also constantly replenishes the hydrogen fuel burnt at the base of the envelope 
by bringing it from its outer parts, where T <^ Tjnax- In general, the results in Fig. show 
how the increase of the T^ax peak value changes the relative abundance distribution. We 
conclude that SPPN nucleosynthesis simulations with nova Tmax-trajectories provide a proper 
qualitative indication of the behavior of their corresponding complete MPPNP simulations. 
However, they can be used only as a diagnostic for nova nucleosynthesis, e.g. in reaction 
rate sensitivity studies. 

The nucleosynthesis i n ONe novae and its sensitivity to reaction rates have been stud- 
ied using post-p rocessing ( llhadis et al.ll2002l ) and full hydrodynamic f lJose fc Hernanzl 120071 : 



Jose et al.ll2010[ ) models. As a result, three reactions have been identified whose rate uncer- 
tainties have the most significant impact on models' predictions: ^^F(p,a)^^0, ^^Al(p,7)^^Si, 
and 30p(p,7)3is. 

The ^^F(p,Q!)^^0 reaction rate is crucial for understanding the most intense 511 keV 



^The Barcelona composition specifies the initial abundances only for isotopes lighter than Ca. 
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7-ray emission that can be observed from novae, which is produced by positron annihila- 
tion associated with the decay of ^^F. This nucleus is destroyed in a nova environment via 
^®F(p,a)^^0. The uncertainty in the rate of this reaction therefore presents a hmit to inter- 
pretation of any futu re observed 7-ray fl ux. This reaction rate has been evaluated recently 
in the compilation of llliadis et al.l (120101 ). where its uncertainty was about a factor of 2 over 
the temperature ran ge characteristic of explosive hydrogen burning in novae. Since then, 
several experiments (IBeer et al.ll2011t lAdekola et al.ll2011t iMountford et al.ll2012l ) have been 
performed to study the nuclear structure of the co mpound nucleus, ^^N e, above the proton 
threshold. As a result, this rate has been updated ( lAdekola et al.ll201ll ). and its uncertainty 
has been reduced to a factor of 1.2 at 0.1 GK. 



The ^^Al(p,7)^^Si reaction rate was evaluated for the flrst time by llliadis et al.l ( 120011 ). 
Understanding this rate is important for an accurate prediction of the yield of ^^Al syn- 
thesized in ONe novae. The isotopic abundance ratio of ^^Al to ^^Al serves as a marker for 
identifying the candidate sources for presolar meteoritic grains. Due to the lack of knowledge 
of the properties of the proton reso nances in ^^Si, the uncertainty in this rate at 0.1 - 0.4 GK 
ranges over 4 orders of magnitude (llliadis et al.ll2002l ). To re duce the unc e rtaint y , the proton 



resona nces in ^^Si have been studied experimentally, e.g. by lMatic et al.l (120101 ). IChen et al. 



(120121 ) , and others, and, as a result, the most recent evaluation of the ^^Al(p,7)^^Si rate by 
ChenI (I2OIOI ) reports an uncertainty of about a factor of 670 over the temperature range of 



novae. 



Finally, the reaction ^''P(p,7)'^^S drives the nuclear activity in ONe novae in the atomic 
mass region above A = 30, and affects the isotopic abundance ratio of ^°Si to ^^Si, which in 
turn helps id entify novae as a potentia l origin for some SiC presolar grains. Its rate was flrst 
estimated by iRauscher fc ThielemannI (I2OOOI) via Hauser-Feschbach statistical calculations. 



(120011 ). The properties of the proton resonances 



and was further discussed in llliadis et al. 

^^S were also largely unknown at the time, and therefore, this rate t oo suffered from a n 



m 



Wrede et al. 


2009; 


Parikh et al. 


2011: 



uncertainty ranging over 4 orders of magnitude for nova temperatures (llliadis et al.l 20021). 
Tti T-ooonf ,rooT.o fino TimoiooT- ofT-mof^vo r^f 31 g ]2as been stu died extensively (iMa et al. 



2007 



has an uncertainty of a factor of ~ 17 at 0.3 GK (jParikh et al.ll201ll ) 



We have calculated the effect of uncertainties in the rates of two of the aforementioned 
nuclear reactions on flnal nova abundances with the NuGrid SPPN code. In order to esti- 
mate the maximum effect, we have used the Two = 12 MK trajectory from the upper panel 
of Fig. US] that has the highest value of T^ax = 344 MK. Figs. IT^ and ITBl show the flnal rela- 
tive abundances of stable isotopes resulting from varying the ^^Al(p, 7)^^81 and ^°P(p,7)'^^S 
reaction rates, respectively, within their (largest) 4 orders of magnitude uncertainties used 
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by 


Iliadis et al. 


(2002) 


Iliadis et al. 


(2OO2I) in 



are in good agreement with those presented by 



abundances are more affected by reducing the rates considered, starting respectively from 
28Si and 30Si. 

The destruction of ^^Al is mainly due to the ^^Al(/3"'")^^Mg and ^^Al(p, 7)^^81 reactions. 
Reducing the ^^Al(p, 7)^^81, the nucleosynthesis path synthesizes less ^^Si and ^^P, therefore 
reducing the flux through this nucleosynthesis channel. On the other hand, ^^Al is accu- 
mulated causing a more efficient decay to ^^Mg. Starting from ^^Mg, the proton capture 
chain ^^Mg(p,7)^^Al(p,7)^''Si(p,7)^^P(/3"'")^*Si causes a higher ^®Si production. However, the 
suppression of the ^^Al(p, 7)^^81 channel makes the proton capture nucleosynthesis overall 
slower, reducing the production in the heavier S-Cl region and increasing the yields of Si and 
P in the final nova ejecta. 

The '^°P(p,7)'^^S reaction is the dominant destruction channel for ^"^P. Therefore, reduc- 
ing its rate makes ^°P a bottleneck for further proton-capture nucleosynthesis, reducing the 
production of heavier species. The isotope ^"^P is accumulated, and also '^°Si via '^°P(/3'^)^°Si. 
As a result, more "^^Si and ^^P are made, and less between S and ^°Ca. 

The differences between the final nova abundances caused by the variations of the reac- 
tion rates are much smaller than those obtained for the slightly different WD's initial central 
temperatures in Fig. [TH Therefore, these nuclear uncertainty studies for novae are particu- 
larly important for detailed comparisons of the isotopic composition in a given atomic mass 
range, e.g. from presolar grain measurements. 

Notice that these tests based on changing a single reaction rate are quite common in 
nuclear astrophysics studies. This needs to be done carefully, since a consistent uncertainty 
study for nuclear reaction rates of the nearby species should also be included. It is not the 
purpose of this section to carry out such a detailed analysis for ^^Al(p, 7)^^81 and ^''P(p,7)'^^S. 
The point that we wanted to stress here is that single-zone trajectories can be used as 
powerful tools for studies to identify which reactions are relevant for the nucleosynthesis 
(i.e., for sensitivity studies), and within which errors reaction rates are needed to not affect 
stellar abundance predictions (i.e., for uncertainty studies). 
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5. Effects Caused by ^He Burning 

In Paper I, we have found, for tlie first time, tliat in tlie extreme case of a very cold 
CO WD, e.g. witli TwD = 7 MK, accreting solar-composition material with a very low rate, 
say M = lO~^^M0/yr, the incomplete pp I chain reactions lead to the in situ synthesis 
of '^He in a slope adjacent to the base of the accreted envelope and that the ignition of 
this ^He triggers convection before the major nova outburst. The ^He burning continues 
at a relatively low temperature, T ^ 30 MK, approximately until its abundance is reduced 
below the solar value, only after that the major nova outburst ensues triggered by the 
reaction ^^C(p,7)^^N. Although this is an interesting variation of the nova scenariojfl because 
the CO enrichment of the accreted envelope is produced by the ^He-driven CBM in this 
case, its observational frequency is expected to be very low (Paper I). Here, we confirm 
this result for ONe novae (Fig. [T7|) . Furthermore, we have found that in the 1.15 Mq ONe 
WD with a relatively high initial central temperature. Two = 15 MK, accreting the solar- 
composition material with the intermediate rate, M = 10~^'^MQ/yr, the interplay between 
the ^He production and destruction in the vicinity of the base of the accreted envelope first 
leads to a shift of T^^x away from the core-envelope interface followed by the formation 
of a thick radiative buffer zone that separates the bottom of the convective envelope from 
the WD surface (Fig. [TSl) . This result is obtained only for the case when an ONe WD 
accretes solar-composition material and it almost disappears in the models with the pre- 
mixed accreted envelopes. Given that the formation of the radiative buffer zone is revealed in 
the computations with the more realistic nova model parameters, such cases can probably be 
observed. We defer a more detailed study of this peculiar case and its possible consequences 
to a future work. In particular, it would be interesting to see if the CBM can cross the buffer 
zone and reach the WD core. 



■^It is interesti ng that '^He was considered as the most hkely isotope to trigger a nova outburst by 
SchatzmanI (|195ll ). 
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6. Comparison With Observed Chemical Compositions of Novae 
6.1. Element Abundances From Optical and Ultraviolet Spectroscopy 

Difficult ies with dete r minin g element abundances from the spectra of nova ejecta are 



discussed by iGehrz et al.l ( 119981 ) . The same authors have compiled the largest set of mass 
fractions of H, He, and heavy elements, such as C, N, O, Ne, in novae from optical and 
ultraviolet spectroscopy that were published in the period from 1978 until 1997. In Fig. [191 
we compare these abundances with the corresponding nucleosynthesis yields predicted by 
ou r l.lSMf;^ nova rnodels from Figs. [11] and [121 We also include the ab undance pre d iction s 



by I Jose fc Hernang (119981 ) . For each of the novae listed in Table 2 of iGehrz et al.l ( 119981 ) , 
we have used all the available abundance measurements, therefore some of the novae are 
represented by up to three points for the same element in our Fig. [TH] that correspond to 
different data sources. We have arbitrarily divided the observed objects into CO (panel A) 
and ONe (panel B) novae using the neon mass fraction A(Ne) = 0.08 for their separating 
boundary value. The panel B shows a good agreement between the observed and predicted 
element abundances for the ONe novae. A comparison for the CO novae shows the presence 
of a large number of observed objects with too low C and O abundances as compared to the 
model predictions. The latter rather represent upper limits for the observed data. We do 
not know how to interpret this discrepancy, other than to assume that the observed material 
has been mixed with the solar-composition material from the WD's companion. We have 
tried to reduce the amount of CO WD's material in the pre- mixed accreted envelope from 
50% to 25% but this does not help much (compare the solid and dashed blue curves in the 
panel A). Note that our model predictions almost coincide with those of the Barcelona group 
in Fig. [H 



6.2. Presolar Grains 



Absorption feature s due to carbide grains have been observed around novae (e.g.. lStarrfield et al 



19971 : iGehrz et al.lll998l ). This opens a possibility to find in the solar system, hidden in car- 
bonaceous pristine meteorites, presolar grains that condensed around nova objects shortly 
before th e Sun formed, and p ossibly carrying the Ne-E component observed in meteoritic 
samples ( IStarrfield et al.l Il997l ) . This scenario is confirmed by the analysis of abundances 
for a small fraction of presolar grains identified so far, which show distinctive isotopic sig- 
natures that may be assoc iated with nova nucleosynthesi s : a minor fraction of t he presolar 



silico n carbides (SiC nova, lAmari et al. 



20071 ). a few graphite grains (e.g.. 



2001 



Amari 



Amari et al. 



2002 



2001 



Jose fc Hernanz 2007; Heck et al. 



Zinner et al.l 120071 ). and oxides (e.g.. 



- 18 - 



Gyngard et al.l l2010t iNittler et al.l l2010l : iGyngard et al.l l201ll ) . Oxide grains should form in 
CO novae, whereas carbide grains around ONe novae, since the ejected material shows the 
C to O abundance ratio s less than unity in the first case, and being C-rich in the last case 
(e.g., lAmari et al.ll2001t also, see our Table However, observational features of carbide 



dust have been observ ed also around CO novae, and therefore this simple distinction cannot 
be made ( lGehrzll2002l ). The distinctive isotopic signatures of nova grains include extremely 
high excesses of ^^C, ^^N, overabundances of ^^Mg and ^^Ne due to the later radiogenic decay 
of the radioa ctive ^^Al and ^^Na (Table [H), and high ^^O/^^O isotop i c ratios, specifically fo r 
oxide grains (IClavton fc Hovlelll976l : lAmari et al.ll2001t lAmarill2002l : I Jose fc Hernanzl 120071 ). 
Nittler fc Hoppe (j2005 ) suggested that, at least, part of the nova SiC grains are instead 
made from material ejected by core-collapse SNe, showing typical nova signatures coupled 
with ^^Si and ^^Ti excesses that cannot be obtained from novae. This is potentially question- 
ing also the origin of other carbide nova grains identified so far, where only C and N isotopic 
ratios were measured. Isotopic abundances from different elements (e.g., C, N, Si, Ca) need 
to be measured in order to disentangle between these two scenarios. On the other hand, the 
fact that carbide grains are observed condensing around novae, makes with no doubt these 
objects a stellar source to keep on the watch in grain abundance analysis. 

In this section, we compare the model results presented in the previous sections for ONe 
novae with presolar SiC grains. Also plotted are the results for CO novae from Paper I. 
We derive a single data point for the isotope abundances from each nova model, using the 
mass average of all zones in the last time step of post-processing nucleosynthesis simulations. 
This is a common procedure used to compare stellar model predictions with isotopic ratios 
measured in presolar grains. We do not distinguish between carbon and oxygen rich layers in 
the expanding nova envelope. This results in an average isotopic composition for each indi- 
vidual nova simulation and does not necessarily represent what presolar grain incorporate. A 
more detailed comparison between grains of possible nova origin and ejected layers with nova 
models would require a dedicated investigation. Figs. [20] and HH compare presolar SiC nova 
grains with our model results. For presolar grains, we used the Washington University (St. 
Louis) database for data extraction (IHynes fc Gyngardll2009l ). The data mostly fall between 
the model and average solar system isotopic ratios. Therefore, present models reproduce the 
data well, if some contamination with solar system material is taken into account. Sample 
contamination either happened on the asteroid parent body or in the lab. Since no major 
and minor element concentrations were measured in the grains, a mixing calculation cannot 
be performed here. 

Equilibrium condensation calculations show that SiC grains are only ex pected to con - 
dense in environments where the C/0 ratio is larger than one (see, e.g., IZinnerl |2003| ). 



Therefore, unless mixing of the ejecta takes place, the chosen approach to mix all shells 
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together, is not completely valid. Fig. [22] shows the evolution of the C/0 elemental ratio in 
the ONe nova models (upper panel) and CO nova models from Paper I (lower panel). To 
make changes on the surface visible and plot all models in one graph, the mass layers were 
normalized from zero to one and inverted, such that the outermost shell plots on the left 
side. Only part of the nova ejecta is C-rich. Without a clear mixing prescription for the 
ejecta - if existing - a definite comparison with presolar grains requires a careful study where 
the analysis is performed grain-by-grain, and considering the isotopic signature in different 
stellar zones. 

In summary, we find a general agreement between our nucleosynthesis calculation trends 
and isotopic abundance measurements from SiC grains with a possible nova origin. The 
comparison however only gives a trend of our nova models since mixing in the ejecta and 
contamination with solar system material have to be further taken into account and dis- 
cussed. A more detailed grain-by-grain analysis is required for a consistent and constraining 
analysis, including analyses of more isotopic constraints where available. 
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7. Conclusion 

We have created the Nova Framework that allows to simulate the accretion of H-rich 
material onto a white dwarf (WD) leading, for a suitable set of initial parameters, to a nova 
outburst, and to post-process its accompanying nucleosynthesis. The Nova Framework com- 
bines the state-of-the-art stellar evolution code MESA and post-processing nucleosynthesis 
tools of NuGrid. It includes a number of CO and ONe WD models with different masses 
and central temperatures (luminosities) that can be used in simulations of nova outbursts. 
The use of the Nova Framework is facilitated by a number of shell scripts that carry out the 
routine job necessary to coordinate the operation of the MESA and NuGrid codes. 

Within the Nova Framework effort, we provide a large set of nucleosynthesis calculations 
for nearly 50 CO and ONe nova models at the solar metallicity [Z = 0.02) so far. To verify 
our calculations, we compare their results with those published in the literature for similar 
nova models. The comparison shows a very good qualitative agreement analogous to the case 
when different nova yields from the literature are compared with each other. Typical features 
of abundance signatures of novae from previous studies and observations are confirmed. For 
instance, the large production factors for ^^C, ^^N, and ^^O in both CO and ONe novae, the 
synthesis of ^Li in CO novae, the overabundance of Ne as a distinctive feature of ONe novae, 
as well as the accumulation of relatively large mass fractions of the radioactive isotopes ^^Na 
and ^^Al in ONe novae (unfortunately, not detected by 7 telescopes yet). 

The nucleosynthesis in novae is studied using the post-processing method. We show that 
such a technique provides abundance predictions consistent with stellar model calculations. 
We study the impact on the final abundances of model parameters relevant for novae, such 
as the mass, composition and initial central temperature of the underlying WD, and the 
mass accretion rate. In particular, decreasing the initial central temperature causes the 
outburst peak temperature to increase, which in turn enhances the production of isotopes 
in the region between Si and Ca. The main effect of the decrease of the accretion rate is 
an increase of the envelope mass Mace accreted before a nova outburst. The higher Mace 
value results in a stronger explosion with a higher peak temperature. For accretion of 
solar-composition material with a very low rate, e.g. M = lO~^^M0/yr, onto a very cold 
WD, e.g. with Two = 7 MK, we find the in situ synthesis of ^He taking place via the 
incomplete pp 1 chain near the bottom of the accreted envelope. The ignition and burning 
of this '^He at a relatively low temperature (T ~ 30 MK) triggers convection in the envelope 
before the major nova outburst, the latter ensuing only after the '^He abundance has been 
reduced below its solar value. We demonstrate that this happens not only in CO novae 
(Paper 1) but also in ONe novae. Moreover, we reveal that the interplay between the ^He 
production and destruction in the solar-composition envelope accreted with an intermediate 
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rate, e.g. M = lO~^°M0/yr, by the 1.15 M0 ONe WD with a relatively high initial central 
temperature, e.g. Two = 15 x 10^ K, leads to the formation of a thick radiative buffer zone 
that separates the bottom of the convective envelope from the WD surface. 

We use three Tmax-trajectories, extracted from full nova models, in single-zone post- 
processing nucleosynthesis computations with the NuGrid SPPN code. We show that the 
results obtained in the one-zone simulations are qualitatively consistent with the yields from 
complete models post-processed with the NuGrid multi-zone code MPPNP, including the 
impact of stellar parameters, like the WD initial central temperature. Therefore, simple 
trajectories can be used as important diagnostic tools for nova nucleosynthesis, e.g., for 
nuclear sensitivity and uncertainty studies. As example, we have studied the impact on nova 
nucleosynthesis of the ^^Al(p, 7)^^81 and '^*^P(p,7)^^S reaction rates. In particular, we find 
that abundance predictions are more affected by reducing those rates. 
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Fig. 1. — Upper panel: Tmax-trajectories extracted from our MESA 1.15 Mq CO nova sim- 
ulations for TwD = 10 MK (log^g -^^wd/-^© = —2.69; for the correspondence between WD's 
initial central temperature and luminosity, Two and L^vdj see Tabled]) and M = 10~^^ Mq/jt 
carried out using the small (nova.net, 33 isotopes from H to ^^Mg coupled by 65 reactions, 
solid red curve) and extended (nova_ext.net, 48 isotopes from H to ^°Si coupled by 120 
reactions, dashed blue curve) nuclear networks. Lower panel: relative (with respect to so- 
lar) abundances (mass fractions) of stable isotopes calculated with the NuGrid multi-zone 
post-processing nucleosynthesis code MPPNP, averaged over the indicated mass range in the 
expanding envelope, in a final model of our 1.15 M0 C O nova simulations with Two = 15 
MK and M = 2 x lO~^°M0/yr (compare with Fig. 1 of I Jose fc Hernanzl (119981 ) and Fig. 5 
in Paper I). 
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Fig. 2. — Evolutionary tracks for pre-mixed CO (solid red curve) and ONe (dashed blue 
curve) nova models computed with the MESA code for the same initial parameters: Mwd = 
1.15M0, TwD = 12 MK, and M = 2 X IQ-^^Mq/ji. 
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Fig. 3. — Initial abundances (solar-scaled mass fractions) of isotopes in the pre- mixed ma- 
terial accreted by our 1.3 Mq ONe WD model (red circles) are compa red with those used in 



similar nova models by t he Barcelona group (Jose et al. 1999 . 200l|) (gree n squares in the 



upper panel) and by both Politano et al. ( 19951 ) and 



Starrfield et all (120091) (blue diamonds 



in the lower panel). In all the cases, it is assumed that the WD accretes a mixture consisting 
of equal amounts of its core and solar-composition materials. Element symbols point to their 
corresponding most abundant isotopes in the Solar system (^^C, ^^N, ^^O, ^"^Ne, ^^Na, ^^Mg, 
27 Al, etc.) 
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Fig. 4. — Tinax-trajectories from two ONe nova simulations with the same values of Two = 12 
MK and M = 2 X 10"^° Mq/jt, but for different WD masses: Mwd = 1-15 Mq (panels A and 
B), and Mwd = 1-3 Mq (panels C and D). The trajectories are compared with one another 
for different nuclear networks and initial compositions (Fig. [3]). The largest network used 
in our simulations is nova_weiss.net. It includes 77 isotopes from H to ^"^Ca coupled by 
442 reactions. If either network or abundance set is not indicated for a curve, it has been 
computed with our standard ONe nova network nova_ext . net or MESA initial abundances, 
respectively. 
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Fig. 5. — Upper panel: the net of 147 isotopes (H, He, Li, Be, and B isotopes are not 
plotted) used in our post-processing nucleosynthesis calculations. Reaction fluxes (arrows) 
correspond to our SPPN simulation for the dashed blue trajectory from upper panel of Fig. [8] 
to a moment when the temperature reaches its peak value of Tmax = 408 MK (Table [1]). 
Middle and lower panels zoom out the isotopes and reaction fluxes for the NeNa-, MgAl-, 
and CNO-cycles. 
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Fig. 6. — Upper panel: Tmax-trajectories from two 1.3 Mq pre-mixed ONe nova models with 
^WD = 20 MK and M = 2 x 10^^^ Mq/jt computed using the same largest nuclear network 
(nova_weiss .net) but different MESA reaction rate data, NACRE and JINA Reaclib. Lower 
panel: solar-scaled final abundances for these two models extracted from their MESA output 
files. Compare the red symbols from the lower panels in this figure and Fig. [TJ the latter 
having been obtained in the post-processing nucleosynthesis computations with NuGrid's 
MPPNP code. 
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Fig. 7. — Tmax-trajectories (the upper panel) and solar-scaled final abundances (the lower 
panel) from two 1.3 Mq ONe nova models, using the largest (the solid red curve and red 
circles) and extended (the dashed blue curve and blue squares, the latter being overlapped 
by the red circles) networks. For these models. Two = 20 MK, and M = 2 x 10~^° Mq/ji. 
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Fig. 8. — Same as in Fig. m but for Two = 7 MK, and M = 10'^^ Mq/ji. For such extreme 
cases leading to the highest values of Tmax? the switching to the nova_weiss.net nuclear 
network is necessary. 
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Fig. 9. — Upper panel: solar-scaled final abundances in the expanding envelope of our 1.3 Mq 
ONe nova model with Two = 12 MK a nd M = 2 x 10_~^° Mp-J yr. The initial abundances for 
these simulations have been taken from 



Politano et al 



, ^ - ( 119951) . Lower panel: the ratios of our 

final abundances to those reported by Starrfield et al. (j2009 ) for a similar nova model with 
MwD = 1.35 Mo, TwD = 12 MK, and M = 1.6 x IQ-^^ Mq/jt (I2005A in their Table 5), 
who used the same initial abundances. 
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Fig. 10. — Same as in Fig. |9l but here we liave used the Barcelona initial abundances and 
compare our MESA /NuGrid simulation res ults with those reported and plotted for a similar 
ONe nova model by I Jose fc Hernang (119981 ) in their Table 3 (the 0Ne6 column) and Fig. 3. 
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Fig. 11. — A comparison of Tmax-trajectories (the upper panel) and final abundances (the 
lower panel) from two of our 1.15 ONe nova simulations with Two = 12 MK and M = 
2 X 10~^^ Mq/jt. The solid red curve and red circles are the results obtained with the 
50% pre-mixed initial abundances in the accreted envelope, while the dashed blue curve and 
blue squares represent the case with the solar-composition accreted material and convective 
boundary mixing (CBM) modeled using the same method and CBM parameter /nova = 0.004 
that we used in the 1.2 Mq CO nova model in Paper I. 
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Fig. 12. — Same as in Fig. [TTl but for the corresponding CO nova model. 
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Fig. 13. — A comparison of Tmax-trajectories (the upper panel) and solar-scaled final abun- 
dances (the lower panel) from three of our 1.3 Mq pre-mixed ONe nova models computed for 
the same accretion rate, M = 2 x 10"^*^ Mq/jt, but different WD initial central temperatures 
(the symbols in the lower panel have the same color coding as their corresponding curves in 
the upper panel). The initial isotope abundances are those of the Barcelona group. 
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Fig. 14. — Lower panel: the final solar-scaled abundances of stable isotopes calculated with 
the SPPN code using the trajectories from the upper panel of Fig. [13] that correspond to 
the three different WD's initial central temperatures: Twd = 12 MK (diamonds), 15 MK 
(squares), and 20 MK (circles). Upper panel: the abundances from the first two sets are 
plotted as ratios with respect to the abundances from the last set. The unstable isotopes 
were allowed to decay. 
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Fig. 15. — Comparison of the final abundances of stable isotopes calculated with the SPPN 
code using the T-wd = 12 MK trajectory from the upper panel of Fig. [inifor the recommended 
set of reaction rates (Xother) with the abundances (Xi) obtained with the rate of the reaction 
^^Al(p,7)^^Si increased (cir cles) and decreased (diamonds) by the factor of 100, like in the 
nuclear sensitivity study of llliadis et al.l (j2002[ ). 
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Fig. 16. — Same as in Fig. [151 but for the reaction ^°P(p,7)^^S. 
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Fig. 17. — A snapshot of profiles of relevant stellar structure parameters in the envelope 
of our 1.3 Mq ONe nova model with T-wd = 7 MK and M = 10^^^ Mq/jt at the moment 
following the ^He ignition that has triggered convection. Like in the CO nova model with 
the cold WD and low accretion rate discussed in Paper I, the low Two and M values favour 
the accumulation of "^He in a slope adjacent to the core-envelope interface followed by its 
ignition at a relatively low T, before the major thermonuclear runaway ensues. 
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Fig. 18. — A snapshot similar to that shown in Fig. [T71 but for our 1.15 Mq ONe nova model 
with TwD = 15 MK and M = 10^^" Mq/jt. The interplay between the ^He production and 
destruction shifts the peak temperature away from the core-envelope interface, located at 
the leftmost steps on the abundance profiles (the middle panel), to a place where the ^He 
burning generates the maximum energy. This leads to the formation of a convective zone 
separated from the interface by a buffer zone. Note that this model and the one in Fig. [17] 
accrete solar-composition material. 
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Fig. 19. — Comparison of the H, He, C, N, O, and Ne mass fract ions in novae from optical 
and ultraviolet spectroscopy (red star symbols, iGehrz et al.lll998l ) with the theoretical ones 
obtained in nova simulations. Blue circles connected by solid blue curves in panels A and B 
correspond to our CO and ONe nova models with 50% pre-mixed accreted envelopes from 
Figs. [12] and [m The dashed blue curve in panel A represents our results for a 25% pre-mixed 
CO nova model, while green squares connected by dot-dashed gree n curves in panels A an d 
B are data for nova models, C05 and 0Ne3, from Tables 4 and 3 of I Jose &: Hernang (119981 ). 
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Fig. 20. — Comparison of the C and N isotopic ratios from our ONe and CO nova models 
with pres olar grain da t a. Bl a ck hnes indicat e the s o lar system va l ues. P r esolar grain data 
are fr om kmari et al.1 fcoOlh : ICao fc Nittleil JiQQtI i: Ifloppe et all fll996h : iNittler fc Hoppe 
(120051 ). Nova models are denoted by their composition (Ne or CO), WD mass (Mq) and ini- 
tial central temperature (MK), followed by a number of models in the evolutionary sequence. 
Their expanding envelopes are mixed as described in the text. 
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Fig. 21. — Comparison of the Si isotopic ratios from ONe and CO nova models with sin- 
gle presolar grain data. For the isotopic ratios, the delta notations are used ((5(ratio) = 
(stellar ratio/solar ratio - l)xlQOO). B l ack li n es indicate the solar system value s . Th e 
presolar grain data are fromlAmari et al.l (|2001[): ICao fc Nittlerl fll997l ) ; iHoppe et all ( l2010l ) ; 



Nittler fc Alexander! (120031 ) : lNittler fc Hoppd (120051 ). Nova models are denoted by their com- 



position (Ne or CO), WD mass (Mq) and initial central temperaturei (MK), followed by a 
number of models in the evolutionary sequence. Their expanding envelopes are mixed as 
described in the text. 
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Fig. 22. — C/0 ratios in expanding envelopes of all nova models in the last time step. The 
C/0 = 1 is indicated by a straight line in each plot. To plot all models on the same mass 
scale and point out differences in the outer layers better, the x-axis shows the mass-shell for 
each model normalized from to 1 and inverted. This means that the outermost shell lies 
on the left side, while the innermost shell is on the right. The figure clearly demonstrates 
that mixing of the whole envelope is not a priori valid, rather carbon rich zones have to be 
mixed and analyzed in order to obtain a valid comparison with presolar SiC grains. 



Table 1. Parameters of Computed Nova Models 



WD 


Mwd[Mq] 


rwD[io' K] 




Afacc[lO-5M0] 


lgiH[Lo] 


T,„ax[10« K] 


Nucl. Network 


Init. Comp. 


CO 


1.0 


12 


-2.55 


5.2 


10.88 


200 


nova 


MESA 


CO 


1.0 


20 


-1.99 


2.1 


9.59 


163 


nova 


MESA 


QQa,h 


1.15 


10 


-2.69 


4.0 


11.84 


257 


nova, ext 


MESA 


CO 


1.15 


12 


-2.50 


2.8 


11.27 


236 


nova 


MESA 


CO" 


1.15 


12 


-2.50 


9.7 


11.28 


243 


nova 


MESA 


CO 


1.15 


15 


-2.25 


1.6 


10.45 


208 


nova 


MESA 


CO 


1.15 


20 


-1.94 


0.94 


9.67 


185 


nova 


MESA 


ONe*" 


1.15 


12 


-2.54 


4.0 


10.24 


262 


nova, ext 


MESA 


ONe*" 


1.15 


12 


-2.54 


4.0 


10.14 


263 


ext, weiss 


MESA 


ONe'^ 


1.15 


12 


-2.54 


15.0 


9.89 


245 


ext 


MESA 


ONe 


1.15 


20 


-2.03 


3.0 


9.88 


245 


ext 


MESA 


ONe 


1.15 


20 


-2.03 


2.0 


9.46 


219 


ext 


Barcelona 


ONe^ 


1.3 


7 


-3.05 


3.0 


10.82 


408 


ext, weiss 


MESA 


ONe 


1.3 


12 


-2.52 


1.7 


10.23 


357 


ext 


Politano 


ONe 


1.3 


12 


-2.52 


1.5 


10.60 


344 


ext 


Barcelona 


ONe 


1.3 


15 


-2.29 


0.95 


9.90 


313 


ext 


Barcelona 


ONe 


1.3 


20 


-1.98 


1.0 


9.86 


318 


ext, weiss 


MESA 


ONe 


1.3 


20 


-1.98 


0.52 


9.26 


267 


ext 


Barcelona 


ONe 


1.3 


20 


-1.98 


0.53 


9.29 


266 


weiss (NACRE) 


Barcelona 


ONe 


1.3 


20 


-1.98 


0.52 


9.37 


265 


weiss (JINA) 


Barcelona 



^This model uses the accretion rate M = 10"" M^/yv, while aU other Af = 2 x 10"^" Mq/jt. 

^For simulations done for two different nuclear networks, e.g. nova_ext.net and nova_weiss.net, only the results for the 
first case arc shown (for the network highlighted with a bold font). 

'^This simulation includes CBM with / = 0.004, like in Paper I, while all other assume that the WD accretes a mixture of 
equal amounts of its core and solar-composition materials. For this case. Mace gives the mass of H-rich envelope that includes 



WD's outer layers penetrated by the CBM. 
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Table 2. Comparison of Selected Abundances in Pre-Mixed Nova Models' 



Abundance Fig. [U COS Fig.\n\ 0Ne3 Fig.^ 0Ne6 



C/0 0.56 0.78 0.36 0.43 2.2 1.2 

X22 6.1 X 10-^ 2.9 X 10-^ 2.6 x lO^* 5.3 x 10"^ 1.1 x 10"^ 6.0 x 10~* 
X2e 2.1 X 10--* 4.7 x 10"^ 3.1 x lO"'* 9.3 x 10"* 8.4 x lO"'^ 7.2 x lO""* 



''The C to O number ratios and mass fractions of ^^Na an d ^^Al are taken f rom o ur models 
shown in the indicated figures and from Tables 3 and 4 of Jose fc Hernanz (199^ for their 
corresponding counterparts that are denoted by C05, 0Ne3, and 0Ne6. 



